A REMARK ON EXACT FORMULAS FOR THE RIESZ ENERGY OF THE 

NTH ROOTS OF UNITY 
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. Abstract. In the paper the author, Hardin and Saff [Bull. Lond. Math. Soc. 41 (2009), 

^v^j ' no. 4] we provided the complete asymptotic expansions of the Riesz s-energy of the JVth roots 

of unity. In this note the exact formulas (valid for all N > 2) are obtained for s is an even 
integer. Explicit representations of the coefficients of the powers of TV are derived by direct 
computations. In the case of positive s a continuous modified energy approximation of the 
Riesz energy is used. Stirling numbers of the first kind, Eulerian numbers and special values of 
partial Bell polynomials play a central role. A number of identities between these quantities are 
[**~- , shown. 

(N 

The Riesz s-energy (real s ^ 0) of an iV-point configuration Xn = {zi, . . . , Zn} in the complex 
Qh| plane (identified with the Euclidean space M 2 ) is given by 

N N 

j=i k=i 1 3 fci 

It can be shown (see [3]) that the iVth roots of unity (or any rotation thereof) maximize (if 
—2 < s < 0) and minimize (if s > 0) the above energy functional over all iV-point configuration 
on the unit circle. We remark that for s < —2 an optimal configuration on the unit circle consists 
of half of the points at one of two opposite locations. (This follows from results of G. Bjorck [3].) 
jy-^ The Riesz s-energy (if s ^ is real) of the iV-th roots of unity can be evaluated by means of 

>/~) ■ N-l , , x -s 

(1) £ S (N) = 2- S NJ2 sin TT ' ^V>2, seC, s^0. 



-I— > 



> 

o 



It has the following general complete asymptotic expansion of C S (N) ([H Thm. I.I]) 

2 p 

(2) £ S (N) = V S N 2 + J2 a n(-s)as-2n)N 1 + s - 2n + O s , p (N- 1+ncs - 2 P), as N -> oo, 

[Zl: > n=0 

provided s is not a positive odd integer. (These exceptional cases are considered in [4l Thm. 1.2].) 

The Riemann zeta function has its trivial zeros at negative even integers. Thus, the sum in ([2]) 
vanishes for all N > 2 if s is a negative even integer or can have at most m terms if s = 2m is a 
positive even integer. 

The case s = —2m and m a positive integer. In [4] Remark 1, Eq. (1.19)] we observed that 

C S (N)=V S N 2 , s = -2,-4,-6,..., 

given in which is correct up to a remainder term of order SjP (N~ 1 ~ s ~ 2p ), where 
p can be any positive integer. The precise result is as follows. 

Proposition 1. Let m be a positive integer and set s = —2m. Then 

m / \ 

£ s (N) = V s N 2 + 2N 2 Y / (-V k [ _J, # = 2,3,4,..., 



fc=i 

N\k 



v m — kj ' 

where N\k denotes that the integer N divides the integer k. The sum above vanishes if N > m. 
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For the sake of completeness we include here the proof. 
Proof. Using the identity (0 Eq. 1.1.9]) 



m— 1 / \ /0 



(sin0) 2m = 2 1 - 2 " 1 ^ (-l)™- fc ( ™* ) cos[2(m - fc)0] + 2- 2m 
in ((T|) we get 

A.- — Q ^ — 1 \ A, — J 

a representation in terms of trigonometric sums satisfying the identity ([8j Eq.s 4.4.4.6-8]) 

N—i 

(3) cos(2irk£/N) = 
Also observe that [4J Eq. (1.8)] can be written as 

(4) Vs = 2 - a r((l - s)/2) = 22m r((l + 2m)/2) = 22m r(m + 1/2) T(l + m) 



N UN divides £ (that is N\£), 
otherwise. 



/7pr(l-s/2) Fr ( 1 + m ) V^Fr(l + m)T(l + m) 

T(2m + 1) _ /2m> 
~ r(m + l)T(m + 1) ~\m, 

by the duplication formula for the gamma function. The result follows. □ 
The case s = 2m and m a positive integer. In [4] Remark 1, Eq. (1.20)] it was observed that 
2 - 

(5) C S {N) = — -Y,an(s)((s~2n)N 1+s - 2n , s = 2m, m = 1, 2, 3, . . . , 

^ ' n=0 

by arguing that in the general asymptotic expansion ([2]) the constant V s given in ((4]) vanishes if s 
is a positive even integer and the Riemann zeta function part £(s — 2n) assumes the value for 
s — 2n is an even negative integer. One should add that (|5|) is true up to a remainder term of 
order O s , p (N~ 1 ~ s ~ 2p ), where p can be any positive integer. In fact, direct computations 
suggest that relation (J5J) holds for all N > 2. The coefficients a n (s) arc defined by the following 
generating function relation and can be expressed explicitly in terms of generalized Bernoulli 
polynomials Bj° (x) ; that is 
(6) 

smirz \ s ^ 2 „ (-l) n B { 2 s J{s/2) 2n 

= 1^ a n{s)z , \z \< 1, s e C, a n (s) = -— 2tt , n > 0. 

ttz I L — ' (2n)\ 

For non negative even integers the Riemann zeta function can be expressed in terms of Bernoulli 
numbers by means of ({2m) = 2 2m - 1 (-l) m - 1 B 2m Tr 2m /(2m)\, m = 0,1,2,.... 

The validity of relation ([5]) for all N > 2 follows by means of a limit process from the corre- 
sponding result for a modified (and thus non-singular) form of the Riesz energy of the iVth roots 
of unity. 

The distance between a point x on the unit circle (|x| = 1) and a point y inside the unit circle 
( r = |y| < 1) is given by 

|x - y| 2 = 1 - 2r(x,y) + r 2 = (1 - r) 2 + 4r [sin(</>/2)] 2 , cos 4> = (x,y). 



This motivates the following modification of the the s-energy of the JV-th roots of unity: 

^ / 2-Kk \ ~ s/2 

(7) X s (A;r):=A^ l-2rcos— +r 2 , < r < 1. 



k=l 
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Observe that C S (N) = lim^i- A4 s (N;r). It should be mentioned that Chu [5] considered sums 
of the form 



n-1 2 



fm(z) = Y^ 



[l-2zcos(27rfc/n) + z 2 ] m ' 

where m and n are positive integers and z is complex. The following recurrence relation is given: 

, / \ fm(z) , z d / fm{z) 

Jm+l\Z) 



1 — z 2 m dz 1 1 — z 2 



Chu's main interest was in the case when z = y— T. 
We make use of the Gauss hyper geometric function 



oo 

(8) 2 F i(°^;*)=E 



(a) « (&) « i*i < i, 



n=0 ^« nl 

where (a) n = r(n + a)/T(a) denotes the Pochhammer symbol. 
Proposition 2. Lei s > 0. Then for integers N > 2 

oo 

X S (7V; r) = N 2 G {s; r)-N(l- r)~ s + 2N 2 £ G^(s; r), < r < 1, 

i/=i 

where the series is uniformly convergent with respect to r on compact subsets of [0, 1) and 
G n (s; r) = (1 - r 2 ) 1 -* 2 F X (1 " s/2. ^ 

For s a positive even integer, that is s = 2m, we can write (cf. [T| Eq. 15.8.7]) 

The hypergeometric function above reduces to a polynomial of degree m — 1 . 
For s = 2, that is m = 1, one simply has G„(2; r) = r™/ (1 — r 2 ) which leads to 

-/V 2 AT 2.V- A ,. v A 2 TV 2N 2 r N 

21 ;r, ~ 1-r 2 (1-r) 2 + 1-r 2 ^ r ~l-r 2 (1 _ r ) 2 + (1 - r 2 ) (1 - r N ) 

l + r N TV 



iV z - 



(1 _ r 2) (1 _ r iV) - {1 _ r y 

Using Mathematica, a limit process r — > 1~ yields that 

w 3 AT 2 fix 2 ir 2 

C 2 {N) = lim M 2 {N:r) = — -- and by ©: £ 2 (iV) = — j —AT 3 - —TV 
r-^-i- 12 12 47r z \ 6 6 

It follows that (J5J) is correct for s = 2 and all TV > 1. 

In the general result we make use of the Stirling numbers of the first kind s(n, k) defined by the 
expansion ([Tj Eq. 26.8.7]) 

n 

(10) ^2 s ( n > k ) xk =x(x-l)---(x-n + l) = (x + l- n) n , 

k=0 

where (a) n is the Pochhammer symbol, and the Eulerian numbers (J£ ) given by (cf. [TJ Eq.s 26.14.5 
and 26.14.6]) 



fc=o x 7 j=o 

' \ /n> 



There holds ( ° > = (™) = 1 and ( n " x ) = 1. By convention we set ( ™ ) = for k > 



n. 



4 



J. S. BRAUCHART 



Proposition 3. Let m and N be positive integers. Then 



f ~\t \ n2 ^r(2m-fc-l) 
M 2m (N;r) = — — 



r(m) ^— J T(m — k)k\ 



k 

N\-q- 1 



-(1 - m) k + 2^2g(k,q;N, 1 - m) (l - r 

g=0 



(1 - r 2 ) 



2 \ fc+l-2m 



7V(l- r )~ 2m , 0<r<l, 



where 

g; a, 6) = V(-l) p s(^P + 95 &)&(P + <7,PK +g , 6(p + ?,p) = V 

— — \ 7 / V £> 



k—q q . , 

p + q\ fp + q-j 



s(n, £; y) = ^ I ) s(n, fc) (y + n - l) fc ~ 



We are interested in the limit as r — > Clearly, .M 2m (iV):= IuHth-i- -M-2m(N\ r) = £2 ra (JV). 
Proposition 4. Let m and iV fee positive integers. Then 

N 3 N 2m 

(12) M 2 (iV) = — --, X 2m (7V)=^/3 2m _,(m)7V 1+2m ^, to>2, 

where the coefficients /3 v (m) (0 < v < 2mJ do wot depend on N. 

Explicit expressions for the coefficients j3 v (m) are given in the Proof of Prop. 2] (see Eq.s (fTTH) ). 
The expansion ([12")) holds for every integer N > 2, whereas the expansion ([5]) (derived from an 
asymptotic result as N — > oo) representing the same quantity holds for sufficiently large N. Since 
the coefhcients of the powers of N in either expansion do not depend on N, by comparing (|19l) 
and ([5]) we obtain for s = 2m and to a positive integer the connection formulas 

2((s — 2n) 

P 2n +i{m)=0 (0 < n < m), ^2m-2n(w) = — t— - — a„(s) (0 < n < to). 

(Zn) b 

The coefficients a n (2m) can be expressed in terms of generalized Bernoulli polynomials and the 
Riemann zeta function for nonnegative even integers can be expressed in terms of Bernoulli num- 
bers (cf. ([H]) and text there). This leads to 

fe - 2 " (m)= (2m -2n)! (2^)1 ' 

For n = one gets the following identity. 

Proposition 5. For m = 1, 2, . . . on has 

(-IT-'B^ ^ " 2 2 +- 2 -r(2m-l-^) ^ 

= V M" 1 tv rnT7 ^ b(v,p)G(2m — p,2m — p,v — p+1), 

where the numbers b(v,p) and G(n,n,q) are given by 

^ = E ( • ) ( k ~ p j ) . £ ( ~ irG(n ' n ' 9) *" = (^) - L 

1. Proofs 

The Gegenbauer polynomials Cn of degree n with index a > may be defined by means of 
the generating function relation 

(13) f>"C( A )(cos0) = - — ^- -x, |z|<l,A>0. 

„=o (l-2zcos0 + z 2 ) 
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They admit a representation through trigonometric functions (cf. Section 22]) as follows 

" (X) (X) 

(14) (cos <j>) = Vfn_V C ° S[( " " m] ' A f °- 

We need the following auxiliary result. 
Lemma 6. Le£ s > 0. Then 

A n (s;N):=Nj:C^(coB^) = N^ £ (s/ gf_ 7 ^ - AT C^d). 
fc=l f=o ^ 

For even and odd n one ftas 

AT|(2£) 



/I f„. ATI _ 9A r2 ( s / 2 ),y-£( s / 2 )„+l+f r(s/ 2 )m/v 

A 2 „ +1 ( S ,iV)-2iV 2^ + -c 2v+1 (i)i\r. 

iV|(2«+l) 

Proo/. By flU) and © we have 

, , m ^ (s/2) t ( S /2) n _ t 2nk " (s/2) e (s/2) n _ e ( , /2) 

fcO v ' fc=l fco v ' 

A|(n-2£) 

Since N\(n — 2£) if and only if N\(2£ ~ n), where 21 — n = — [n — 2(n — £)}, we can use symmetry 
to write 

(s/2) e (s/2) 2v _ e ^ (s/2) e (s/2) 2 „_ e (s/2)„(s/2)„ 

ei(2v-£)\ f^ Q l\{2v-l)\ v\v\ 

N\(2v-2t) N\2(v~l) 

2 ^ (s/2Us/2) 2u+1 _ e » 

^ a(2u + l-£)\ ^ £\{2v+l-t)\ ■ 

N\(2v+1-%1) JV|(2w+l-2i) 

The results for even and odd n follow by reversing the order of terms on the right-hand sides above 
and substituting the expressions into the result for general n. □ 

Proof of Proposition [H Using (1131) in ([7]) and Lemma [5J we obtain 

M s (N;r) = ±r-Nj:C^(co S 2 -^)=^±^ £ ^^ - N £^^(l). 

n=0 k=l n=0 1=0 >' n=0 

N\{n-1l) 

The last series has the closed form (1 — r)~ s (cf. (fT3")0 . One has 

V 2, v ^fZgMfZ 2 ^ = y (fZgWg), 2, | 2 y 2, y ( S /2)^( 5 /2), +£ 

S ^ *M ^ (*-*)!(./ + 

jV|(2i/-2£) JV|(2«) 



y (f/Wg), 2 - (s/2) v _ t (s/2) u+l 



1=1 
N\(2t) 



V 

(2t)l ' Zii V 2^+1 



2 F 1 (^/ 2 ;r 2 )+2 £ i^fV&J'/W + '/lS 



(i 
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In the last step © was used. By the same token 



E^ +1 E 



N\(1v+l-2l) 



(s/2) t (s/2) 2v+1 _ e 
£\(2v + l-i)l 



(f/^Wi 2£+i F fs/2,2£ + l + s/2 2 

2^ + 2 



2 ^ (2^+l)! r 



2-T 1 



AT|(2£+1) 



Putting everything together we obtain 



M S (N; r) = N 2 G (s; r)-N (1 - r)-'+2N* f] G„(s; r), G n (s; r):=^^r B 2 F X ( s / 2 > n + s / 2 : r 2 



n=l 
Win 



71+1 



The series J2^Li G u n(s; r) converges uniformly with respect to r on compact subsets of [0,1), 
which can be seen from the integral representation (cf. [TJ Eq. 15.6.1]) and estimate 



G n (s;r) =r" 



r(n + s /2) r 1 W 2 - 1 (l - *) B — /a 



[T(s/2)] 2 r(n + 1 - fl /2) 7 (1 - r2i)" +s / 2 



T(n + 5 /2) „ 
< c s (r)— : tttt r 



c s {r) 



dt 



t s/2-l 



■dt, 



T(n + l-s/2) ' "* v '' [r(s/2)] 2 y (l-r 2 t) 6 
valid for n > no > s/2, where c s (r) can be uniformly bounded on compact sets in [0, 1). Hence 

r(n + s/2) 



0<J2 G »N(s;r)<C+ J2 G n {s-r)<C + c s {r) ^ 



^ r(n + l - s/2) 

n— no 



r™ for some C > 



i/=l n—riQ 

and the right-most series above converges uniformly on compact sets in [0, 1) for n > s/2 > 0. 
Application of the last linear transformations in [TJ Eq.s 15.8.1] yields 

G n (s; r) = (1 - r 2 ) ^ 2 F X (1 " * A "+ 1 - s/2. ^ 

For s a positive even integer, that is s = 2m, we can write (cf. [TJ Eq. 15.8.7]) 

-r" a Pl( 1 -^tJ~ m ;l-r 9 y 



G»(2m; r) = (l - r 2 ) 1 " 2 " 1 (m) " ■ ' ^ -»'■" + ] " '"■ i - 



where the ratios can be simplified further. This shows ©• □ 
Proof of Proposition The series expansion of the hypergeometric polynomial in (J9j) yields 

V^ r M £^ (2-2m) fe /c! 

Since 



k=0 



(1 - m) fe _ fm - 1\ I (2m - 2\ _ r(m) T(2m - 1 - fc) 
(2 - 2m) fe ~ V k )/\ k ) ~ T{m-k)T{2m- 1) 

and r(2m — 1) = 2 2 ™ l ~ 2 r(m — 1/2) r(m)/y/ir (duplication formula for the gamma function) we 
obtain 



G „ (2m;r) . _±_,,| E^a^o. (1 _ ^ 



2 m 



By Proposition [2] and the last representation we get 

T(2m-k- 1) 



M 2m (N;r) = 



N 2 rn-1 



E 



r(m) ^— ^ r(m — fc)fc! 



£ (HJV + l-m) fc rl"l 



;i-r 2 ) fc+1 - 2m -7v(i_ r )" 2 ™. 



By Lemma [15] the infinite series reduces to a rational function in r . That is 



M2m(N-r) 



N 2 



m—1 

E 



Tpm- k - 1) 
r(m) ^-^ r(m — fc)fc! 



-(l-m) fe + 2^. 9 (fc,(?; AT,l-ro)(l-r 

9=0 



iVN-9- 1 



(I-, 2 ) 



2 \ fc+l-2m 



JV(l-r) 



-2m 
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where g(k, q; N, b) is defined in Lemma [T5l □ 

Proof of Proposition [4} For the Proof of Proposition 2] we need some preparations. We define 
the following two functions 

/l_ r AT\-8 ^ 

/«W = (t^t) ■ W = ^ = E ^ 

v 7 ' 1=0 

and get a series expansion of f q (r) at r = 1. For that we use Faa di Bruno's differentiation formula 

' K=l ^ ' 7 fe=l 

where fc = fci + ■ ■ ■ + k n in the first sum and this sum is extended over all partition of n, that is 
integers ki, ■ ■ ■ ,k n > such that fej + 2k 2 + ■ • • + nk n = n. The polynomials B nt k(xi, X2, • • ■ ) in 
the second sum are the (partial) Bell polynomials, explicitly given by 

^ B nAxuX2 ,...)= £ t^t(tt) (#) < 

fei + fc2H — =k 
fci+2/c2H — — n 

which satisfy the generating function relation (see, for example, [6]) 
(16) 



°° fin i / °° y.m \ 

^B n ^ lll2 , l3 „..)- = - Yi^mU ' 

n— \m— 1 / 



(For the polynomial B ni fc only the variables x\, X2, ■ • ■ , x n ^k+i will be effective.) We record that 
(17a) B n , k (axi,ax 2 ,...) = a k B nt k(xi,x 2 , . . a^O, 

(17b) B n:k (a 1 x 1 ,a 2 x 2 ,-.-) = a n B n ^(x 1 ,x 2 , . . .), a ^ 0. 

Lemma 7. Let n be an integer > 0. Then 

w n + 1 \« + V 

Proof. One has 

jy n (n) (N-l ) W JV-l JV-1 



/l (»)(r) = <|L^-} = ^£r4 =E^} W = E / ( / - 1 )-('- n )^ 

=0 J «=o 



1 -r 



with the understanding that a sum with upper index is smaller than the lower index is empty and 
therefore set to zero. If r = 1, then 

AT-l N—n N-n , . , , 

^)(i)=En^i)---(^-^)-E^+"- i )---^+ i )^Ew»= ( 1 } ; 1 ■ 

f=n £=1 fcl 

The last step follows by complete induction in N. □ 
Lemma 8. Le£ q be a positive real number and N a positive integer. Then for m = 1,2,... 

/,(r) =N- q + J2 ll ^T 1 ~ !)" + K m (r), < r < 1, 

n— 1 

where the coefficients fqiX) can be written as 

= D-i)^)^-* B «^ 1! ( 2 ) ' 21 (3) ' 3 ' ^ 

and £/ie remainder can be estimated by 



/-. \m-f-i m-t-A / A7\ / A7\ / A7 

n m (r)\ < V^pryr E («)*/«+*(»•) B m+ i, fe (i! ( 2 ) , 21 ( 3 ) , si ' - * 
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Proof. Let m > 1. By Taylor's theorem 
f q (r) 



m f (")/ 1 x ,(m+l), x 

, , , ^AJl) (r _i)» + 72m(r ), ^(r) = ^ T ^(r-ir +1 

' n— 

for some p with r < p < 1. By Faa di Bruno's differentiation formula 

O) n 



I / 1 — " \ I . . 

4 n) (r-) - | (t^7 j } = E Bn^'M, ft"(r), . . . ) 

n 

= E(- 1 )"(9) fe Wr))"*"* B„, fc (fc'(r), fc"(r), . . . ). 



fc=i 



The representation of fq (1) follows from the fact h(l) — N and Lemma [7] 
The remainder lZ m (r) is estimated next. Observe that the function h^ u \r) is positive and 
strictly monotonically increasing on (0, 1) for each v > 0. Thus, by the definition of Bell polyno- 
mials (fT5l> one gets 



m+l 



f { q m+1) (p)\ < £ {q) k {h{p))-«- k \B m+l , k {h'{r),h"{r),...)\ 

k=l 
m+l 

< E («)* (Mr))-* - * B ro+1)fe (^(l), h"{l), ...). 



k=l 

The estimate follows. □ 
Let ( aH hafc ) denote the multinomial coefficient. 

The specific Bell polynomials appearing in Lemma [8] can be represented as follows. 
Lemma 9. Let n, k and N be positive integers. Then 

v 7 v 7 v 7 ni — 1 rik—l 

niH hn fc — n 

k-\-n 



n 



k\(n + k)\ L 



J2(-l) n - i+k H e (n,k) N e 



n -\- k 



where the coefficients 

n n 

H ti n, »-£...£ ( Bi + . . +JE- E w». + moi ■ ■ ■ w« + 

ni=l n*=l N 7 li=l t k =i 
niH \-n k =n h^ fc =f 

vanish for I = 0, 1, . . . , fc — 1 and do not depend on N . 

Proof. By the generating function relation for partial Bell polynomials (|16[) 

f-w..©.K?).K^)....,S4(|,U l H t 
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The right-hand side above evaluates as 



1 

jfc! 



(l+tf - 1- Nt 



-i k 



1 



N 

E 



M=2 

N-l N-l 

~ E ' ' ' E 

1/1 = 1 v k = \ 

k(N-l) 



-i k 



N 
vi + 1 



1 

ft! 



N-l 

E 



N 
v+l 



-i k 



N 
v k + 1 



v—k 



E { k\ E " ■ E 



JV-l iV-l 



"1=1 ft = i 

l/l H hl/fe=l/ 



N 
v x + 1 



ii/iH 



N 
Vk + 1 



t v 



This shows that the infinite series at the left-hand side above is, in fact, a finite series and com- 
parison of coefficients yields the first identity of our result 



B n , fc (l! 



2! 



,3! 



) 



kl 



N-l N-l 

E-E 

ni = l n k =l 
m-\ hnfc=n 



N 
ni + 1 



N 
n k + l 



Let [x] n denote the falling factorial [x] n = x(x — 1) • • • {x — n + 1). Taking into account that 
(„+]} = \N\ v+ i_/(v + 1)! vanishes if v + 1 > TV for the positive integer N, we may write 



B„, fc (l! 



2! 



,3! 



fc! ir + K + 

riiH \-nu=n 



By the definition of Stirling numbers of the first kind and using that n\ H + n k = n one has 



(m + 1)! K + l)! 



=n 



i 



(n„ + 1) 



1=0 



ni+l "k + l 

E-E 



s(m + 1,4) • • -s(«fe + 1,4) ^i + ... + 4 



<i=0 

E- 



4=0 

rii + l rifc + 1 

E-E 

4=0 4=0 
(, 4 + ...+4=£ 



(m + !)!••• (rife + 1)! 

s(ni + 1,4) • • • «(^fc + 1,4) 



(m + !)!■■■ (rife + 1)! 



When expanding the left-hand side above in powers of iV we infer that the lowest power appearing 
at the right-hand side above is N k . That is 



ni+l n k +l 

E-E 

4=0 4=0 
£i+---+e k =i 



s(ni + 1,4) • • • s(nfc + 1,4) 
(m + !)!■■■ K + 1)! 



= for I = 0,1,..., A;- 1. 



Putting everything together we get the identity 



B„, fc (l! 
and the representation 



2! 



,3! 



e=k 



n roi+1 n k +l 



W (n,fc)= E ••• E E-E 

m = l n fc =l 4=1 4=1 
niH hn fc =n 4H h4=£ 



s(ni + 1,4) ■ ■ -s(n k + 1,4) 

(m + !)!■■ -K + l)! ' 



in 
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It is well-known that (—1)™ s(n, fc) > for n > k. Thus 

k fc 



J] s(n„ + 1,4) - |s(n„ + 1,4)1 = 

i/=i i/=i 

and from n\ + • • • + n k = n and 4 + • • • + 4 = £ it follows that 

Bm(H ( 2 ) , 2! ( 1 ) , 3! ( 1 ) , . . • ) = g E (-1)-^ *)| ^. 

By definition (— l)"~ fc s(n, fc) counts the permutations of 1, 2, ... ,n with exactly k cycles. With 
this combinatorial interpretation in mind we can write 

where the coefficients He(n, k) (which do not depend on N) are given in the lemma. □ 

Special values of the numbers Hg(n, k) are as follows. 
Lemma 10. Let n (> 1) and k (> 1 ) be integers. Then 

H n+k (n,n)=2- n (2n)\[ 



n 

n+l n+l 



v — • \ — > / n + k \ kl(n + fc)! „ ,11 
H n+k (n,k)= E""E = 1 Mn.o.-- 

rvH i^ n fc — n+A; 

Proof. By Lemma [9] 

n n ✓ _ \ ni^l n n -\-l 

Mm)=E-E L lB8 + ; „ +1 E-E + 1,4)1 -|s(n n + 1,4)1 

niH hn„=n €xH h^n=n+& 

It follows that ni = ---= ri„ = 1. Hence, the multinomial coefficient evaluates as 2~™(2n)! and 

1 1 

H n+k (n,n) = 2- n (2n)\ E E |s(2,4 + 1)| • • • |s(2,4 + 1)| • 

<iH h£n=k 

Since s(2, 1) = —I and s(2,2) = 1, the sum counts the number of n-tuples (4,..., 4) with 
precisely fc of the entries being 1 which is equal to (?). The result follows. 
Similarly 

n n / , \ ni+l n k +l 

H„ +fe (n,fc)=E-E L +1 n ll n k + l) E-E Kni + l,4)|--|s(n fe + l,4)|. 

niH h«fc=n |-£ fc= ri+fc 

Here, the inner sum reduces to a single term with 4 = "i + 1, • ■ • , 4 = "fc + 1 which evaluates 
to 1. Thus 

71 71 / , s 71 + 1 71 + 1 , , 

H n+k { n ,k) = y ... y ( n + k )=y...y n + k 

TllH h«fc=« 7liH h7lfc=7l+fc 

On the other hand, by Lemma [9] 

n! 

H n+k (n, fc) 



fc!(n + fc)! 

is the coefficient of the highest power of N in 
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From 

v - („ + 1) = ^Ti NU+1 ^ + ° {N ^ as N 00 

and the identities (fl"7|) we obtain 
It follows that 

n 1 11 

^-—^^(n, k) = B n , fe (-, -,...), n > k > 1. 

Lemma allows us to recast the representation of fq(l) given in Lemma |5] 
Corollary 11. Let q, n and N be positive integer. Then 

e=o k=i ^ 

where the numbers Hi + k{n,k) are given in Lemma\^ 

Proof. From Lemmas [8] and [9] 

. n J k-\-n 

n 1 n 

fc=l ' V £=0 

= E(-!)' E TuS^r g ^( n ' fc ) 

Z_. ^ k\{n + k)\ 

The result follows. 

Lemma 12. Le£ q be a positive integer. Then one has the generating function relation 

E(-l) n G(n,n l(7 )x"- (!— -1. 



n=l 



Proof. By Corollary HT1 and Lemma [TU] 

(-l)"G(n, n, ff ) = E fc7 (n + fc) ! g ^^ k ^n\ D" 1 )*^ 1/3, . . . ). 
Application of the identities (cf. Faa di Bruno's differentiation formula) 

g (/w) = E Efc=1 B -* (oi>a2 -- ) »", /(*) = EV- ^) = EV> 

z — ' 7i! z — ' n! ^ — ' n\ 

n—1 n—1 n—1 

to the functions 

/(^ErS^^^ ^HE^^/ = (i + ,r-i 

n=l v ; n=l 

yields 



E(-l) n GKn,< Z ).T"= 3 (/(.T)) = 



Ti=l 



- 1. 
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With these preparations we are ready to give the following proof. 

Proof of Proposition [^} The formula of M.2m(N] r) in Prop.[3]is rearranged w.r.t. powers of (1— r). 
That is 



i i / »t \ , at /i \-2rn N' 2 \ -» r(2m — fc — 1) Nfc+l-2rn /, xfc+l-2m 

M 2ro fiV; r) + (1 - r) = ~fR g r(m - k)k\ (1 " m) * (1 + r) (1 " r) 

+ 2 FR g g r(m- fe ) fc ! g ^ « ^ 1 " W) (T^r J (1 + r)* 1 (1 - r)* « 
Note that k + 1 — 2to < 0. Using Taylor expansion (cf. Lemma |S1 since f q (r) = (1 + r)~ q if N = 2) 

(m + 1)! 



(l + r)- q = y;^2-»-*(l-r) T, + W TO (r), £ m (r) < ? %±i- (1 - r) m+1 
' 77,' (m + 1 II 

n=0 



one gets (AT is assumed to be fixed) 



, , , s-2m N' 2 r ^T(2m-k-l), 
M 2m (N; r) + N (1 - r) = g - m) fc 

x 2 "U fc (2m-l -fc) n 2fc _ n+1 _ 2m _ fc+n+ i_ 2m + _ r)) 

^— ; 77! 

+ 2 rgy g g r r™ m Z%i ] ^ g; ^ 1 - TO) (t^t) 9 

2m+g-fc . _ _ . 
x £ ^ i 2 fc-»+l-2 TO (1 _ r)n +k- q -2m + _ r)) 

n=0 n ' 
i\t2 m— 12m— 1 — fc „, „ , n 

- F ( m ) 2^ r(m-fc)fc!n! 1 )k ( ' 

K > k=0 n=0 v ; 



rfm) ' ' ^-^ rfm — fc)fc!n! ' ' \ 1-r 

v / fc=0 g=0 n=0 

+ 0((l-r)), asr->l" 
By Lemma [5] 



-1 2m-l-fc 



A^ 2m (7V; r) + N (1 - r)" 2 ™ = £ " jf" ^ 2W ."^ 1} (1 - m) fc 2*+— (1 - r)"+ ! ' " 1 ■*»' 

1 (to) z ^ *-^L I (to — k)k\n\ 



fc=0 n=0 



— EE E ( It ] m)^^^ (1 - r)<* *-<-*» 



rfm) r— ' ' rfm — fc)fc!77! 

v ' fc=0 q=0 n=0 v ' 

2 m-1 fc 2m+q-k2m+q-k-n , , , / i \„ <■ (f) /-, s 

lo N r(n + 2m-fc-l (-1J / g+ i(l) 

v ' fc=0 g=0 n=0 i/=l v 7 

x 2 fe+1 - 2m -" (1 - r y+n+k- q -2m + _ r))] as r ^ i-. 
The power (1 — r) -2 " 1 appears only in the triple sum when q = k and 77 = 0. Its coefficient is 
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By the formulas in Prop. [3] and the properties of the Stirling and Eulerian numbers (cf. (|23[1 ) 
g(k, k;N,l-m) = s(fc, k; 1 - m)b(k, 0) N k = s(k, k) N k ^ / k \ = s(k, k)k\N k = k\N k . 

3=0 ^ 3 ' 

Thus, by Eq. 4.2.3.10], 

c_ 2ro = N2 2 ~ 2m £ lf m ~ ( k - l h k = AT2— £ f 2W - 2 " M 2* = 7V 
^— J 1 (to) 1 (to — k) \ to — 1 J 

and, therefore, the term N(l — r)~ 2m could be removed from both sides of the identity above. 
First, we reorder the terms 

— — 

+ 2 ffr E E " rf (2TO ^ fc l +g T + 1 V ^ * 1 " m)AT— 2«— (1 - r)" 
I (m) j^_ _ Q 1 (rn — k)k\(2m + q — k — n)\ 

N 2 ^k2rn^- k n m2m _ 1 _ k)+q _ n + 1) (-1)^(1) 

+ 2 7r-^Z^Z^ 2^~fv Zil^ — T 1 ^9{k,q; N,l - to) f 

1 to ' ^-^ ' I to — fc fc (2m + (7 — fe — n) ^ 

' fe=0 9 =0 n=0 i/=l v ' v y ' 

x 2 n- q -l-2(2m-l-k) ^ _ r ^-r, + _ ^ 

A r2 



V V r(2(m + fc)-n) /-, _ •> 2 «-2(m+fc) _ y 



r ( m ) ^ ^ r ( TO ~ k ) k K m + k-n) 

+ 2— W V V r(2(2m-l-fc) + g -n + l) (-1) (1) 

r(m)^2- 2- 2- r(m-fc)fc!(2TO + q-fc-n)! al ' J (n-v)! 

v ' fc=0 q=Q n=Q v=0 y ' y * ' v ' 

x 2 »-,-l-2(2m-l-*) ( X _ r )-" + 0((1 _ r )) ; as r -> 1~. 

Since lim,.^.!- M.2m{N; r) — A42m(N) and M.2m{N) exists and is finite for all positive integers to 
and N with N > 2 (and iV(l — r)~ 2m can be removed from both sides of the identities above) it 
follows that the cumulative coefficient of a negative power of (1 — r) is zero and we are left with 

M*m(N) = * g ^ + ^, (1 ~ "»)™-i-*2- a(m+fc) 
r(m) ^ r(m - k)k\(m + fc)! m 1 fc 

+ 2 rrrEE ^f (2m ^ r : fc) + +g+ ^ g (^ « * , i - m)^-« 2— ) 

1 (to) ^— ^ 1 (to — fc)fc!(2m + q — fc)! 

iV 2 "^^ 2 'y- fc r(2(2TO-l-fc ) + g -n + l) , (-ir/a(D r H- 3 (,-i- t l 
+ r(TO)^^ ^ r(TO-fc)fc!(2TO + g-fc-n)! ffl ,9 ' ' j n! 

Substituting the expressions for g{k 7 q;N, 1 - to) (cf. Prop.© and /^(l) (cf. Cor-EJ we get 

fATl- N% V T(2(TO + fc)) 2(m+fc) 
• M2m(7V) " "fR ^ r(m-fc)fc!(TO + fc)! (1 _ m >™-i-* 2 
(18) k ~° 

' m—l k k — q m — 1 k 2m+q—kk — q n 

fc=0 g=0p=0 fc=0 <?=0 n=0 p=0 «=0 
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where the primary coefficients are given by 

2 T(2(2to - 1 - k) + q + 1) (-l) p s(k,p + g; 1 - m)b(p + q,p) 



X(k,q,p;m): 



'r(m) r(m - k)k\(2m + q-k)\ 2 9+ 1 + 2 ( 2 ™- 1 - fc ) 



2 r(2(2m -l-fc)+g-rc+l) (-l)r + "+'s(fc,p + g; 1 - m)6(p + g,p) 
q,p, n, i; m):= T{m) T(m _ k)kl(2m + q _ k _ n) < 2 ,+i +2(2m -i- fc )-„ G ("' ^ + 1) 

which depend on the secondary coefficients 

s{n,£;y) = E s(n, fc) (y + n - l) fe_£ , &(p + <7,.p) = E J ' J ■ /v ,; ' 

fc=£ ^ ' j=0 

and G(n, £, g) is given in Cor. QT] (cf. Lemma [S])- 

We reorder the terms in (IT8|) w.r.t. powers of N. The first part is 



m — 1 fc fc— 9 m — 1 k k—p m—lm—lk—p 

EEE «>n m ) NP = E E E m ) NP = E E E "0^- 

fc=0 (?=0 p=0 fc=0 p=0 9=0 p=0 k=p 9=0 

The second part is 

m — 1 fe 2m+q— kk— q n m — 1 k q 2m — q n 

EE E EE^^.^^i^^EEE E E^^-W'^;™)^ 

k=0 q=0 n=0 p=0 t=0 fc=0 3=0 p=0 n=0 £=0 

m — 1 fc 9 2m—q2m—q m — 1 fe fe 2m q 2m q 

= EEEE E X(k,k-q,p,nJ;m)N^ = EEE E E *(M-9,P,M;m)JV<+P 

fc=0 g=0 p=0 f=0 n=£ fc=0 p=0 q-p £=0 ri=l 

m — 1 m — 1 fe 2m — 9 2m — q m—1 m — 1 2m— 9 m — 1 2m — 9 

-EEEE E^^-W^-^l^^EE E E E X(k,k- qi p,n,l;m)N^ 

p=0 fc=p 9=p f=0 n=£ P=0 q=p 1=0 k=q n=i 

m—1 771 m—1 m—1 2m — 9 m—1 2m— p 2m — £ m—1 2m — 9 

EEEEE+E E E E E 

p=0 £=0 9=P fc=9 n=i p=0 £=m+l 9=P fc'=9 n=l 



X(k,k- q,p,n,£; m) N e+P 
X[k,k — q,p,n,v — p\ m) N v 



m—1 v m—1 m—1 2m — 9 2 m— 1 m — 1 m—1 m — 1 2m — 9 

EEE E E + E E EE E 

i/—0 p—0 q—p k—q n=u—p U=m p—u — m q—p k—q ri—v—i 
2m v— m— 1 2m— u+p m— 1 3m— q 
+ E E EEE X{k,k-q,p,n, V -p-m)N». 

u—m+l p=0 9=P fc=9 n=i>—p 

Putting everything together we arrive at 

t\t2 m—1 rfof _i_ h\\ m—lm—lk—p 



iV 



m — 1 1/ m — 1 m—1 2m — q 2m — 1 m — 1 m — 1 m — 1 2m — q 

EEE E E + E E EE E 

i/=0 p=0 9=p fc=q n=v— p z^=m p=z^ — m 9=p fc=9 n=z^— p 



-X"(A;, fc — g,p, n, v — p;m) N" 



2m «— m—1 2m— w+p m—1 2m— 9 

+ iV E E EEE X(fc,A;-g,p J n,i/-p;m)JV*'. 

i/__m+l p—0 q—p k—q n—u—p 
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Direct computations (with the help of Mathematica) give the expected results (cf. (JSJ)) 

N 3 N 

M 2 (N) = , 

V ; 12 12' 

MJN) = N ) + N [ -N N [ +N 1 

y ' 8 V 4 8 J V 8 720 / V 16 72 / V 720 16 

= 1 N N. 

720 720 720 

In general, one has for m > 2 



2m 

l+v 



M 2m {N) = Y J Pu{m)N 1 - 
where for v = 0, 1, . . . , m — 1 and v ^ 1 

m— 1 fc— w y m — 1 m — 1 1m— q 

(19a) /3„(m):=^ ^X(fc,(?, i/;m) + ^^^ ^ X(k,k - q,p,n,v - p;m), 

k=v ij=0 p=0 g=p fc=g n=v—p 

for v = m and m^l 

m-lm-lm-l 2m — q 

(19b) /8 m (m):= ^ ^ ^ ^ X(k,k - q,p,n,m - p;m), 

p—0 q—p k—q n—m—p 

for !/ = to + 1, . . . , 2m — 1 (and m > 2) 

m— 1 ra-lm-1 2m— q 
I3 v {m):= ^2 E E X! x (k,k-q,p,n,u-p;m) 

p=U—m q—p k—q n— u— p 

19c 

i;-m-12m-i;+pm-l 2m— q 

+ E X! X! X! x ( k ,k-q,P,n,v-p;m), 

p=0 q=p k=q n=u—p 

for ^ = 2m (and to > 2) 

m — 1 m— 1 

(19d) /?2m(TO)-= X(k,k — p,p,2m — p,2m — p\m) 

p— fc— p 

and for = 1 (and to > 2) 

" r(m) 4- T(m - k)k\(m + k)\ [ m )™-i-k* 
(19e) k ~° 

v ' ni — lk — 1 1 m-1 m-1 2m— q 

k=l 9=0 p=0 <j=p fc=<j n=l— p 

The result follows. □ 

Next the coefficients /3„(m) are studied in more detail. 

Proof of Proposition^ By the definition of the numbers X{k, q,p, n, t; to) given in the Proof of 
Prop. [4] we get 

m —1 m— 1 

fi2m{m) = X X X(fc, k - p,p,2m - p, 2m - p; m) 

p—0 k—p 

m-1 m-1 2 2+fc-2m rrt _ -| _ jL\ 

= E E (-1)* rfmWrfm-W S(fc ' fc; 1 - "»WMG(2m -P,2m-p,k-p + l). 

p=0 k=p \ J ■ \ ) 
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By the definitions of s(n, k; x) (Lemma IT5]) and b(n,£) (Lemma [H]) we have 
s(fc, k; 1 - m) = s(fc, fc) = 1, 6(fc,p) = ^ ^ * ^ ^ ~ J 
and the numbers G(n,n,q) are defined by the generating function relation (Lemma II 2[) 

oo / X -]\ "~ 

£(-l)»G(n,n,?)x" = (^^) -1. 

n=l V 21 / 

The result follows after interchanging the two summation signs. □ 
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Appendix A. Auxiliary results 
Lemma 13. Let n be a nonnegative integer. Then 

[ £ j s{n,K) (y + n- i' k ' 

£=0 k=i 

Proof. Using (fpQ|) and the Binomial theorem, one gets 



+ y) n = E s ( n ' e > y) xi > s (™> e ' y) := E ( £ ) s ( n ' k )(y + n - !) A 

0— n i — d \ / 



( x + V) n = ({x + y + n-l) + (l- n)) n = ^ s ( n , k)(x + y + n-l) k 



k=0 

n k 



-.k-t 



= EE s ( n > fc )uK (y + n ~ 1 ) k 

k=0£=0 ^ ' 

The result follows after reordering the sum. □ 
Note that 

n 

(20) s(n, 0; y) = ^ s(n, k) (y + n - l) k = (y + n - 1 + 1 - n)„ = (y) n . 

k=0 

Alternatively, one can use the identities ([?] is the unsigned Stirling number of the first kind) 



oo / \ n 
(« + fo )n = E • )«-A b h « = E[ 



j=0 XJ 7 fe=0 

to obtain the representation 



x 



s(ri,i;j/) = E 





"A;" 


) 





(y) n - 



We need the poly logarithm function defined by 

OO 

(21) M*) = Eip l z l <L 

fe=X 

In particular, one has for ^ = — n and n a positive integer the relations (cf. [7J, [S]) 

oo n — 1 / v 

(22) L_n(«) = g^ fc - (1 _ E Q** 1 "^ 1*1 < L 



A REMARK ON EXACT FORMULAS FOR THE RIESZ ENERGY OF THE NTH ROOTS OF UNITY 17 

Lemma 14. Let n be a positive integer. Then 



L-nW- {1 _ z)n+1 

Proof. By (|2"21 and the Binomial theorem 



3=0 



(1-z) 



n-l 

— J2( n -) [i-(i-*r ! 



n+1 



3=0 



(1-Z) 



n— 1 n— J 

n+l E 51 

j=0 1=0 



n\ n — ] 



3 \ I 



(-i)'(i-*r 



Reordering of the sum yields 

1 



ra-1 



L_„(z) 



(i-zr 



3=0 



6=1 



E 

3=0 



n\/n-j 
3 \ *■ 



(-i)'(i -*r 



The sum over the Eulerian numbers gives n! (cf. [TJ Eq. 26.14.10]). The result follows. 



□ 



Note that (using ( q ) = for q > 1) 



(23) 



3=0 



n-l 

E 

3=0 



= r>l 



Lemma 15. Let k be a non-negative integer and a, b complex numbers. Then 

oo k 

(M« + b) h zW = -(b) k +2"£g(k,q;a,b) (1 - zp" 1 , 

V— — oo <j— 

where (s(n,£;y) is defined in Lemma 1731 and b(n,£) is defined in Lemma [T^P 

fc— g 

ff (Af, g; a, 6):=^(-l) p S (fc,p + g; 6)6(p + 9,pK +9 . 

p=0 

Proof. Let /(z) denote the inifinte series above. For fc = one gets 

2z 1 + z 



z < 1, 



/(*)= X! z H =l + 2$>" = l+ r 



z 1 — z 



1/— — OC 



Let fc > 1. By Lemma fT3l one has 

oo k oo k 

f(z) = (1^+2 J> + 6) k z" = (&) fe + 2^ S (fe ) p;&KX>V / = (6), + 2^s(fc,p;6)aPL_ p (z). 

z^— 1 p=0 y—1 p—0 

By Lemma [HI and (fT0|) (and using Eq.s ([20} and (J23J)) one gets 

/(z) = (6), + 2 s(fc, 0; 6) L (z) + 2 £ s(fc,p; b)a ^ ! + l^t^M £ ~ ^ = (6)fc + 2(5)fe _^ 



P =i 



(1-z) 



p 



2 ^ s(fc,p; b)pla? (1 - z)^ 1 + 2 ]T ]T s(fc, p; 6K(-l) £ 6(p, £) (1 - z)-^ 1 ^ 
p=i p=i i=i 
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Simplification of the first three terms and summing up over equal powers of (1 — z) in the double 
sum yields 



/(*) 



-(b) k + 2j2s(k,p;b)plaP(l-z) 



-p-i 



fc-i 

=E 

9=0 



' fc 

E 

:p=l+q 



s(k,p;b)a"(-iy-%{p,p-q) 



(1-z) 



-q-l 



fc-1 



-(b) k + 2s(k,k;b)k\a k (1 - z)' 
-(b) k + 2s(k,k;b)k\a k (1 - z)^ 1 



fc-1 
^E 

9=0 
fc-1 

^E 

9=0 



s(fc,g;%!a« + £ s(k,p; b)aP(-lf-%(p,p - q) 

p=l+9 



(1-z) 



-9-1 



' fc 

E 

.P=9 



s(k,p;b)aP(-lf-%(p,p-q) 



(1-z) 



-9-1 



fc 

2E 

9=0 
fc 

2E 

9=0 



' fc 

E 

'k — q 

E 

.P=0 



S (fe,p;6)a J, (-l) p -«6(p^-g) 



(1-z) 



-9-1 



(-lfs(k,p + q;b)b(p + q,p)a^' 1 



(1-z) 



-9-1 



where we used that (") = for n > 1 and [[25]) , 

References 



□ 



[1] Digital Library of Mathematical Functions. 2010-05-07. National Institute of Standards and Technology from 
http: / /dlmf.nist.gov/ 

[2] M. Abramowitz and I. A. Stcgun, editors. Handbook of mathematical functions with formulas, graphs, and 
mathematical tables. Dover Publications Inc., New York, 1992. Reprint of the 1972 edition. 

[3] G. Bjorck. Distributions of positive mass, which maximize a certain generalized energy integral. Ark. Mat., 
3:255-269, 1956. 

[4] J. S. Brauchart, D. P. Hardin, and E. B. Saff. The Riesz energy of the A^h roots of unity: an asymptotic 

expansion for large N . Bull. Lond. Math. Soc, 41(4):621-633, 2009. 
[5] W.-C. Chu. Summations on trigonometric functions. Appl. Math. Comput., 141(1):161-176, 2003. Advanced 

special functions and related topics in differential equations (Melfi, 2001). 
[6] M. Mihoubi. Partial bell polynomials and inverse relations. J. Integer Seq., 13(4): Article 10.4.5, 8, 2010. 
[7] J. Milnor. On polylogarithms, Hurwitz zeta functions, and the Kubert identities. Enseign. Math. (2), 29(3- 

4):281-322, 1983. 

[8] A. P. Prudnikov, Y. A. Brychkov, and O. I. Marichev. Integrals and series. Vol. 1. Gordon & Breach Science 
Publishers, New York, 1986. Elementary functions, Translated from the Russian and with a preface by N. M. 
Queen. 

[91 E. W. Weisstein. Dirichlet L-Series. From MathWorld-A Wolfram Web Resource. 



http://mathworld.wolfram.com/Polylogarithm.html Nov. 29 2010. 



J. S. Brauchart: School of Mathematics and Statistics, University of New South Wales, Sydney, 
NSW, 2052, Australia 

E-mail address: j . brauchartSunsw . edu . au 



